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Abstract. Take a linear ordinary differential operator 0(2;) = X)i=i 
with polynomial coefficients and set r = maxi—i k{degQi{z) — i). If 
satisfies the conditions: i) r > and ii) deg (3^.(2) = k + r we call it a non- 
degenerate higher Lame operator. Following the classical examples of E. Heine 
and T. Stieltjes we initiated in |6] the study of the following multiparameter 
spectral problem: for each positive integer n find polynomials V{z) of degree 
at most r such that the equation: 

^{z)S{z) + V{z)S{z) =0 

has a polynomial solution S{z) of degree n. We have shown that under some 
mild non-degeneracy assumptions on T there exist exactly spectral poly- 

nomials Vn^i{z) of degree r and their corresponding eigenpolynomials Sn,i{z) 
of degree n. Localization results of [6] provide the existence of abundance of 
converging as n — > 00 sequences of normalized spectral polynomials {Vn,i^ (z)} 
where Vn^i^{z) is the monic polynomial proportional to Vn,i„{z). Below we 
calculate for any such converging sequence {Vn^i^{z)} the asymptotic root- 
counting measure of the corresponding family {Sn,i„ (z)} of eigenpolynomials. 
We also conjecture that the sequence of sets of all normalized spectral polyno- 
mials {V„^i{z)} having eigenpolynomials S{z) of degree n converges as n — > 00 
to the standard measure in the space of monic polynomials of degree r which 
depends only on the leading coefficient Qki^)- 



1. Introduction 

A generalized Lame equation, see e.g. ptT is the second order differential equation 
given by: 

Q2i^)^ + Qii^^)^ + V{z)S = Q, (1.1) 

where Q2iz) is a complex polynomial of degree I and Qi{z) is a complex polynom 
of degree at most / — 1. If we fix the polynomials Q2{z) and Qi{z) then the classical 
Heine-Stieltjes spectral problem, [5], [IS] asks to determine for any given positive 
integer n all possible polynomials V{z) such that p.ip has a polynomial solution 
S{z) of degree n. Such V{z) are referred to as Van Vleck polynomials and their 
corresponding polynomials S{z) are called Stieltjes or Heine-Stieltjes polynomials. 

The next fundamental proposition announced in [5] was the starting point of the 
Heine-Stieltjes theory. (Notice that throughout this paper we count polynomials 
V{z) individually and polynomials S{z) projectively, i.e. up to a non-vanishing 
constant factor.) 

Theorem 1 (Heine). If the polynomials Q2{z) and Qi{z) are algebraically inde- 
pendent, i.e. their coefficients do not satisfy an algebraic equation with integer 
coefficients then for any integer n > there exists exactly ("+'^2\ polynomials 
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V{z) of degree exactly (l — 2) such that the equation has and unique (up to a 

non-vanishing constant factor) polynomial solution S{z) of degree exactly n. 

In [5] we generalize problem p.ip to high order operators as follows. Consider 
an arbitrary linear ordinary differential operator 

c.(z) = ^Q.(z)^, (1.2) 

1=1 

with polynomial coefficients. The number r ~ maxi^i fc(deg(5i(z) — i) is called 

the Fuchs index of d{z). The operator is called a higher Lame operator if its 
Fuchs index r is nonnegative. In the case of the vanishing Fuchs index d{z) is 
usually called exactly solvable in the physics literature. The operator d{z) is called 
non-degenerate if degQk{z) = fc + r. 

Given a higher Lame operator d{z) consider the problem of finding for each 
positive integer n all polynomials V{z) of degree at most r such that the equation 

t)(z)5(z) + F(z)5(z) =0 (1.3) 

has a polynomial solution S{z) of degree n. 

Following the classical terminology we call (|1.3p (higher) Heine-Stieltjes spectral 
problem^ polynomial V{z) a (higher) Van Vleck polynomial, and the correspond- 
ing polynomial S{z) a (higher) Stieltjes polynomial. To move further we need to 
formulate two main results of [B] , see Corollary 1 and Theorem 9 there. 

Proposition 1. For any non-degenerate higher Lame operator d{z) there exist and 
finitely many (up to a scalar multiple) Stieltjes polynomials S{z) of any sufficiently 
large degree. Their corresponding Van Vleck polynomials V{z) all have degree ex- 
actly r. 

The next localization result is of a special importance to us. 

Theorem 2. For any non-degenerate higher Lame operator d{z) and any e > 
there exists a positive integer N^: such that the zeros of all its Van Vleck polynomials 
V{z) possessing a Heine-Stieltjes polynomial S{z) of degree n > as well as 
all zeros of these Heine-Stieltjes polynomials belong to Couvq^. Here Couvq^, is 
the convex hull of all zeros of the leading coefficient Qk{z) and Couvq^ is its e- 
neighborhood in the usual Euclidean distance on C. 

Remark 1 . Proposition \T\ together with Theorem [2] show that the zeros of all Van 
Vleck polynomials having Stieltjes polynomials of degree greater or equal to a suf- 
ficiently large positive integer n are confined in some fixed e-neighborhood Couvq^ 

of ConvQ^. Therefore, there exist plenty of converging subsequences {Vn,i^{z)} 
of normalized Van Vleck polynomials. Here Vn^i^{z) is a monic polynomial pro- 
portional to Vn.i„ {z) and Ki,i„ (z) is some Van Vleck polynomial having a Stieltjes 
polynomial of degree n. 

It seems natural to pose the following two questions. 

Problem 1. What happens with the set {Vn.i{z)} of normalized Van Vleck polyno- 
mials having a Stieltjes polynomial of degree exactly n when n — > oo ? 

Problem 2. What happens with the subsequence {Sn,i„{z)} of Stieltjes polynomials 
whose corresponding sequence {Vn^i^{z)} of normalized Van Vleck polynomials has 
a limit? 

At the present moment we do not have even a conjectural answer to Problem 1. 
Some initial steps in this direction can be found in [T?'. Problem 2 however has a 
satisfactory answer reported below. 
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Definition 1. Given a probability measure /i supported on some subset of C we 
define its Cauchy transform as 

Obviously, is analytic in the complement to the support of /i. 

Definition 2. Given a polynomial P{z) of degree m we define its root-counting 
measure as the finite probability measure given by 
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m 



where j runs over the index set of the set of all distinct zeros {zj} of P{z), 5{z — Zj) 
is the usual Dirac delta function concentrated at zj and kj is the multiplicity of the 
zero Zj of P{z). 

Assume now that a subsequence {Vn,i^ {z)} of normalized Van Vleck polynomials 
of an operator converges as n — > oo to some monic polynomial V{z) of degree 
r. 

Our first result is a far reaching generalization of the Main Theorem of [5] to- 
gether with the main result of [12j describing the asymptotic root distribution of 
Heine-Stieltjes polynomials under the asumptions that one picks a sequence with 
(asymptotically) the same portion of roots in each of the intervals (aj,ai+i). The 
latter process corresponds to the choice of a sequence of Van Vlecks converging (up 
to a scalar factor) to some limiting polynomial. 

Theorem 3. For any non- degenerate higher Lame operator d{z) of order k take any 
subsequence {Vn^i„{z)} of its normalized Van Vleck polynomials converging to some 
monic polynomial V{z). Then the sequence of the root-counting measures of 

the corresponding Stieltjes polynomials {Sn,j„{z)} weakly converges to a probability 
measure fi.g y whose Cauchy transform C^^ yiz) satisfies almost everywhere in C the 
equation 



'^"-<^> = OTP- 

Figures 1 and 2 illustrate the above Theorem [31 
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Figure 1. The union of the zeros of 40 linear Van Vlecks for the operator 
£s{Q{z)S{z)) + V{z)S{z) = with Q{z) = (^2 + l){z - 3/ - 2)(z + 2/ - 3) 
whose Stieltjes polynomials are of degree 39. 

Remark 2. Notice that by Theorem[2]the support of fj,^ y should lie within Convq^ 
for any V{z) which appears as the limit of normalized Van Vleck polynomials. 
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Figure 2. Zeros of 40 different Stieltjes polynomials of degree 39 for the equation 
£j{Q{z)S{z)) + V{z)S{z) = with Q{z) ^ {z^ + l)(z - 3/ - 2)(z + 21- 3). 



Explanation to Figure 2. The smaller dots are the 39 zeros of S{z)] 4 average 
size dots are the zeros of Q{z) and the single large dot is the (only) zero of the 
corresponding V{z). For most of the pictures one observes the typical structure of 
a curvilinear forest with vertices of degrees 1 and 3 only formed by the roots of S{z) 
which connects the roots of Q{z) and that of V{z). At the same time pictures 3-5 
in the second row show the cases of connected support of the corresponding root 
counting measure. 

Our next result describes some properties of a probability measure in case when 
it satisfies the assumptions of Theorem [31 

Theorem 4. Consider a rational function R{z) = where V{z) and Q{z) are 
monic polynomials with degQ{z) — degV^(z) — k > 2. Assume that the exists 
a compactly supported probability measure ^ on C whose Cauchy transform C^{z) 
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satisfies almost everywhere (wrt the standard Lebesgue measure onC) the equation 

Cl{z) ^ (1.5) 



Then such a measure fi is unique (for a given function R{z) ) and its support is 
a curvilinear forest with leaves at the roots ofV{z) and/or Q{z). This support is 

straightened out in the local canonical coordinate w{z) — J^^ \f^^^^' 

The next question seems very natural in view of the latter theorem but we do 
not have even a good guess about its possible answer except for some trivial cases. 

Problem 3. Which rational functions R{z) admit a probability measure as in The- 
orem 

Remark 3. In fact we will first prove Theorem |4] and the use it to settle Theorem [3l 

Acknowledgements. The authors are sincerely grateful to H. Rullgard for im- 
portant discussions of the above topic and to Prof. Andrei Marti'nez-Finkelshtein 
who brought to our attention the crucial notion of quadratic differential with closed 
trajectories. 

2. On probability measure with Cauchy transform whose power is a 

rational function 

This section is devoted to the proof of Theorem ID Below we discuss a number of 
properties of a compactly supported probability measure fj, whose Cauchy transform 
C^(z) satisfies almost everywhere in C the equation 

where V{z) and Q{z) are monic polynomials with no common factors. 

It is clear that C^{z) can have a pole of order at most 1, which means that Q{z) 
can only have roots of multiplicity < fc if (|2.ip is satisfied. Suppose now that Q{z) 
has roots of multiplicity exactly k, say Q{z) ~ {z — ai)'' ■ ... ■ {z — a„i)'^Q{z), where 
Q{z) only has roots of multiplicity < k. Taking the z-derivative in distributional 
sense yields 

V{z)^ 





h '^'KQMJ (z-a^)-...-{z-a,^)dz[Q(,)) ' 

where 6aj is a point mass located at z = a-,- . Proposition [2] below states that the 

(v( )\^^^ 
j must be supported on a finite union of ana- 
lytic curves. If and C~ denotes the one-sided limits of the Cauchy transform 
as we approach the support of /u, then the Sokhotsky-Plemelj formula tells us that 

(C^ (s)—C~ (s)) 

ft = p{s)ds where p{s) = — ^-^ 1^, with ly a properly oriented normal of the 

support, and ds is the usual arc length. Clearly p vanishes only at the roots of 
V{z). From this it immediately follows that the point masses 6a can not lie in the 
closure of the support of p,, since the total mass of fi is assumed to be finite. 



Our goal in this section is to prove that the support of a positive measure /i whose 
Cauchy transform satisfies (|2.ip is a (curvilinear) forrest. By the above remarks we 
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can restrict our attention to measures whose Cauchy transform satisfies (|2.ip with 
Q{z) only having roots of multiplicity < k. The reader will easily check the validity 
of all the arguments given below if we add a finite number of isolated point masses 
to the measure fi. 



Proposition 2. // /i is a probability measure satisfying i2.1]) almost everywhere 
then the support &^ of is the union of finitely many smooth curve segments, and 
each of these curves is mapped to a straight line by the (locally defined) mapping 



Remark 4. Function ^'(z) is often referred to as a canonical coordinate. Proof of 
this Proposition repeats more or less literally that of Lemma 4 in [S] and is included 
for the sake of completeness. We need the following technical statement. 

Lemma 1 (Corollary 2 of [Sj). For a finite set yl C C, a convex domain U and 
a measurable function x '■ U A the claim that ^ > (as a distribution) is 
equivalent to the existence of real numbers Ca, a € A such that xi^) — o, almost 
everywhere in Ga where 

Ga^{z(^ U; ca + ^(az) > Cf, + 3fi(6z); V6 G A}. 

In other words, any (local) domain where x attains a constant value is (locally) an 
angle given by linear inequalities. 

Remark 5 . Proof of Lemma [1] is based on Corollary 1 of [5] claiming that for a 
convex domain [/ C C, a finite set A G C and a subharmonic function v defined in 
U such that 2^ e A almost everywhere one has that v coincides with the maximum 
of a number of linear (non-homogeneous) functions and is, therefore, convex. 

Proof of Proposition^^ Let us first prove it in a neighborhood of zq G 6^ which is 
neither a zero nor a pole of ^(fj- Choose a branch B{z) of ^ in some simply 
connected neighborhood of zq and define ^' as some concrete primitive function 
of B{z) in ri. Let U be some small convex neighborhood of ^'(zo) such that \l/ maps 
some neighborhood of zq bijectively on U. By (|2.ip we can write 

C^(z) = x(*W)SW 

for z € ^~^{U), where x takes values in the set of fc-th roots of unity. Using the 
variable w — "ifiz) one gets 

dc, _dx{^iz))^,_^,_^^,^fdx\ M 



where denotes the puUback of distributions defined in U by the map 5". Since 
the measure /i is positive one gets the relation ^ > which should be interpreted 
in the distributional sense. By the above mentioned Corollary 2 [/ is the union of 
sets Ga whose boundaries are finite vmions of line segments, such that x is constant 
in each Ga- Therefore, (3^n^'^^(C/) '5^^ (^supp is the union of finitely 
many curve segments which are mapped to straight lines by ^E*. 

If zq is a zero or a pole of q|^| , we take a disk D centered at zq and not containing 

any other zeros or poles of gj^j .If 7 is any ray emanating from zq we can define 

single- valued branches of \J'^^ ^nd ^ in D\^. Notice that is continuous up 
to Zq- Let U be any half disc centered at ^'(zo) and contained in "^{D \ 7). Then 
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by the first part of the proof has the required properties in 'J ^{D). Varying 7 
and U, we see that the same holds in a fuU neighborhood of zq. □ 



Lemma 2. Suppose that the Cauchy transform C f^{z) of fj, satisfies i2.1]) and u{z) 
is the logarithmic potential of ^. If is a (locally defined) inverse of ^ which 

is a primitive to \J^7 then the function u o vj/^i coincides with the maximum of 

some number of linear functions and, in particular, is convex (in its domain). 

Proof. Since u is subharmonic we need to check that its derivative belongs to a finite 
set. We prove our lemma in a neighborhood of any zq which is neither a zero of 

V{z) nor a zero of Q{z). Choose a branch B{z) of \J^^ in some simply connected 
neighborhood Vl of zq and define ^'(z) as some concrete primitive function of B{z) 
in n. Let U be some small convex neighborhood of "^{zq) such that Vl/(z) maps 
some neighborhood of zq bijectively on U . We want to show that uo vl/~i is convex 
in U . By (|2.ip we can write 

C^{z) = x{^{z))B{z) 

for z E ^~^([/), where x takes values in the set of fc-th roots of unity. Since 
Cp(z) — 2^ we have using the variable w — '^{z) 

Therefore by the above mentioned Corollary 1 the locally defined function uovj/^i 
is piecewise linear and convex. 

□ 

Corollary 1. // an endpoint of any curve segment in the support &^ of fi is a 
hanging vertex (i.e. not shared by any other such segment) then this endpoint is 
either a zero or a pole of , see Figure 2. 

Proof. If this were false, then take a point p which is a hanging vertex but not a 
zero of either V{z) or Q{z). The Cauchy transform is supposed to satisfy an 
algebraic equation whose branching points are exactly the zeroes of V{z) and Q{z). 
In particular it has no monodromy around p. This implies that the limits of as 
we approach from both sides close to p are the same, which in turn implies that 
p lies off the support of /i. □ 

Proof of Theorem U starts with a series of additional observations. Notice that 
if we fix a branch of \J^^j^ locally near some point which is neither its root or its 
pole and consider a multi-valued canonical coordinate 

globally i.e. take its full analytic extension then ^'(z) will be well-defined and 
univalent only on the universal covering of C \ {Z{V) U ^(Q)) where Z{V) (resp. 
Z((3)) is the set of all roots of V{z) (resp. Q{z)). But due to the existence of a 
measure /i we can choose an almost global representative of ^'(z) on C substantially 
reducing its multi-valuedness. Namely, let f2 be the complement to the support of 
/I, i.e. n^C\ef,. Define *+(z) in n by 

vl'+(z)= / \ogiz-OdKC)- 
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Obviously, is a part of the whole multi- valued function ^'(z) since for any 

z e n one has = Ic = ^t^i^) "^^^^^ ^i^i^) satisfies (|2?T|) . Since ^ is 

never simply-connected the function ^'+(2:) is still multi- valued, namely, going once 
around some connected component K of in the clockwise direction one increases 
'i'+{z) by 27ri x IJ,{K), where IJ,{K) > is the mass of the measure /z concentrated 
on K. Nevertheless the real part 3i{^+{z)) is a well-defined single- valued function 
in Q coinciding with the logarithmic potential u{z) of /i, namely 

u{z) = U{^+{z))^ [ \og\z-C\dfi{0. 

Consider now the family $ of curves in Q defined by the condition: 

= const 

which is obviously independent of the choice of the branch of 4'+(z). One can easily 
show that the gradient grad ^(2:) of u{z) coincides with C^{z). (Here C^(z) is the 
usual complex conjugate of C^{z).) Thus the family <i> consists of the integral curves 
of the vector field Cfj,{z) which is well-defined and non- vanishing in Q.. Moving along 
the trajectories of C^{z) in positive time we increase the value of u{z). Finally, for 
sufficiently large \z\ one has 

Cu{z) ~ and u{z) ~ log|2;|. 
z 

The next statement is very crucial for the proof of Theorem [H 

Lemma 3. 1) Any trajectory ofCf^{z) tends when t — s- +00 either to 00 in CP^ or 
to a root of V{z). 

2) Any trajectory of —Cfj,{z) tends to &^ or, in other words, any trajectory ofCp,{z) 
'starts ' on 6^. 

Proof. According to lemma [5] locally near any point p € 6^ the logarithmic po- 
tential u is given as the maximum of a finite number of linear functions if one 
uses the canonical coordinate. This shows that the gradient flow points away from 
the support of fj, everywhere, except possibly at a root of V{z) where the gradient 
vanishes. At infinity we have C^{z) « |, which easily implies that the point at 
infinity is a sink for the flow defined by C^(z). Denote by the Riemann sphere 
with an e-neighbourhood of {00} U 6^ removed. To obtain the proof we need to 
show that a trajectory cannot stay in a subset of for any e > when t +00. 
By construction any such does not contain singular points of C^(z). Thus, we 
only need to rule out the occurrence of closed trajectories and recurrence in K^. 
But both of these possibilities are indeed forbidden by the fact that C^t(z) is the 
gradient field of a function without singularities in for any e > 0. □ 

Definition 3. A trajectory of C^(2;) ending at a root of V{z) is called exceptional. 

Lemma 4. There exists at most a finite number of trajectories ofC^{z) ending at 
a given root of V . 

Proof. Wlog assume that a root of V{z) lies at the origin and has a multiplicity p. 
Notice that 
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We now have y ■^[fj = zP/^{ap + Op+iz + ...)'^/^ . Choosing a single-valued branch 
of (op + flp+iz + ...)^/'°, say bo + biz + and taking ^ = 2(69 + biz + ...^/p, which 
can be chosen single- valued in a neighbourhood of the origin, we get \J^!^ = S.^^'^ 
which is valid in a neighbourhood of the origin minus the slit defined by the support 
of our measure /i. Integration yields 




Q{t) \ k )^ 



Since there is only a finite number of rays in the ^-plane entering the origin along 
which 5(^(p+*=)/'=) = 0, our claim follows. □ 



Remark 6. The proof of the previous lemma copies some of local studies of zeroes of 
quadratic differentials, see [TB], Ch.2. It is indeed possible to interpret the support 
of /J, as a geodesic of a higher order differential, which in turn allows us to state 
more explicitley exactly how many exceptional trajectories enter a hanging vertex 
of the support of /i. However we postpone detailed study of this topic, see [9j. 

Definition 4. We denote the union of and all exceptional trajectories by 
and call it the extended support of fi. 

The following statement is very essential. 

Lemma 5. The extended support is topologically a tree. 



Proof. Indeed, by lemmas [3] and d] the flow on CP^ \ defined by the gradient 
vector field (z) (which is non- vanishing there) contracts the whole domain to the 
point at infinity. The result follows. □ 



Let be a measure whose Cauchy transform satisfies (|2.1|) and denote U,* = 
C \ T^. Define now the following specialization of the function '^+{z) to Q.*: 

■^++{z)= [ log(z-C)rfM(C) 
Jc 

where z G Jl*. Note that ^'-(-^(z) is multi- valued only up to addition of multiples 
of 27ri. 

Lemma 6. The function \E'++(z) determines a mapping of fi* onto a domain 
H = {w\'Siw > h{^w)} where h is a picewise- continuous function. Moreover, 
^'^^ : H ^ n* is a single-valued function, see Fig. 3. 

Proof. We have already noticed that is defined up to a multiple of 2-Ki 

in Q.* and that ^'^+(2) = C^{z). If 7 is a small curve segment of supp{^) and 
if t/ is a small one sided neighbourhood of 7, then Proposition [2] states that 7 is 
mapped to a straight line segment by ^'++ which is not horizontal and that U is 
mapped to the right of this straight segment. The latter follows if we observe that 
in the proof of Proposition [2] we have x = 1 in C/ and 5R(x) < 1 everywhere. If 
7 is a part of an exceptional trajectory then it is mapped by 4*++ to a straight 
horizontal line segment since exceptional trajectories are level curves of 
Continuing ^ ++ around we obtain a broken piecewise linear curve of the form 
{5ft(w) = bounding a domain H of the form {lk(w) > h{^{w))} where 

/i is a piecewise continuous function. It is clear that ^++ maps onto H with 
boundary to boundary. The function ^{z) = exp(— ^'++(z)) is single-valued on il* 
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and maps ft* U {00} to D = : log|C| < —h{—arg()}, does not vanish in il* 
and has a simple zero at infinity. It foUows that "0 is bijective on fl* U {00} D 
and hence '9^\{w) = ^^^{e^^) is single- valued. This concludes the proof of the 
lemma. □ 




Figure 3. Extended support of a positive measure associated with a rational 
function i^z-a){z-b){z-c)(z-d) ^^^'^ image under the mapping '^^^{z). 

Explanation to Figure 3. The left picture shows the support of a measure /i con- 
sisting of 2 components - one connecting three poles a, 6, c and the other connecting 
the only zero a with the remaining pole d. The exceptional trajectory connects the 
point O on the first component with the zero a so that the extended support 
forms a tree. The map '^++{z) wraps out the complement C \ onto the strip 
bounded by two horizontal lines and the piecewise linear curve shown on the right 
picture. This map has a period 2tt\/—1. We go around the extended support in 
the anti-clockwise direction starting from the pole c. The piecewise linear curve on 
the right picture shows the images of the special points we meet while doing this: 
c, M, O, a, d, a, O, b, M, a, M and returning back to c but with the imaginary part 
increased by 27r. Each time we meet one of the poles or the point M we have to 
change our direction by ±120°. 

We can finally prove the uniqueness of the required measure. Indeed, assume 
that there are two different probability measures ^1 and /i2 whose Cauchy trans- 
forms solve ()2.ip almost everywhere and let 1*1(2;) and U2{z) be their logarithmic 

potentials. Notice that there is only one branch of \J'^^ which has ^ as its asymp- 

totics near 00 in CP^. Therefore the Cauchy transforms and logarithmic potentials 
of /ii and ^2 have to coincide in some neighborhood of 00. We will show that 
U2{z) > ui{z) in Ql = C \ T^^i and ui{z) > U2{z) in Q2 = C \ Tp,2- Indeed, 
Mi(*++(w)) = 'Si{w) for all w e Hi and U2(^'++(w)) = 3?(w) for aU sufficiently 
large ^(w). On the other hand, U2 o is piecewise linear and convex on any ray 
3(u') = const in Hi. Therefore, U2(^'++(if)) > ^(w) for all w G Hi. Changing 
place of ui{z) and U2(z) we get the second inequality. Since ui{z) and U2{z) are 
continuous in the whole C they should coincide. But the measures fii and ^2 are 
obtained as Aui{z) and Au2{z), therefore they coincide as well. □ 
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We finish this section with some information about connected components of 

Lemma 7. Each connected component of 6^ contains either some number of roots 
of Q{z) and equally many roots of V{z) or it contains k + j roots of Q{z) and j 
roots ofV{z) for some j — 0, ...,r. (We count roots with multiplicities here.) 

Proof. Going around such a component should result in a trivial monodromy of the 
branches of \J^^- Notice that these branches are cyclicly ordered and a small 
clockwise oriented loop around a root of V{z) results in the cyclic clockwise shift 
by ^ while a small clockwise oriented loop around a root of Q{z) results in the 
cyclic counterclockwise shift hy □ 

3. Proving Theorem [S] 

Our scheme follows roughly the scheme suggested in [5] . We need to prove under 
its assumptions the sequence {/Xra,i„} of root-counting measures of the Stieltjes 
polynomials {Sn.i„iz)} converges weakly to a probability measure f-xiiz) v '^hose 
Cauchy transform C^^^^^ satisfies almost everywhere in C the equation p.4p . 

To simplify the notation we denote by {Sn{z)} the chosen sequence {Sn,i„{z)} of 
Stieltjes polynomials whose normalized Van Vleck polynomials {Vn,i^{z)} converge 
to V{z) and let {fin} denote the sequence of its root-counting measures. Also let 
fin'' be the root measure of the ith derivative Sn\z). Assume now that NN is a 
subsequence of natural numbers such that 

/iW- hm 

exists for all j = 0,1, ...,k. The next lemma shows that the Cauchy transform of 
ft = p,^^") satisfies the required algebraic equation. 

Lemma 8. The measures fi^ are all equal and the Cauchy trans form C (z) of their 
common limit satisfies the equation (|1.4p for almost every z. 

Proof. We have 

Sk'+'\z) fdpi^ 
{n~j)si'\z)~' J z-C 

with convergence in Lj^^, and by passing to a subsequence again we can assume 
that we have pointwise convergence almost everywhere. From the relation 

i)iz)Sn{z) + Vniz)Sn{z)^0 

it follows that 

Qk{z)Sk^\z) ^ Qi{z) ( ^ Sli+'\z) \ 

n...(n-k + \)Sn{z) n...{n-k + l) {n - l)...{n - k + 1) [j}^ {n ~ j)si'\z) ) ' 

One can immediately check that ~ Ti(n-i)"i;n-fc+i) ~^ while the sum in the 

right-hand side converges pointwise to almost everywhere in C due to presence 
of the factors (n — l)...{n — fc + 1) in the denominators. Thus, for almost all z G C 
one has 

Sl'^\z) V{z) 

n{n - l)...{n -k + l)S'„(z) ^ Qkiz) 
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when n ^ oo and n e NN. 
then one has 



(^) 



,(0) 



(^) 



If u^^'>{z) denotes the logarithmic potential of 



1 



lim — log 

n^oo n 



1 



n{n — l)...{n — k + l)S'„(z) 



lim - log|y(z)| 

n^oo Tl 



loglQfe(z) 



0. 



On the other hand we have that u^'^^z) > m*^^)(z) > ... > u^'^^z), see Lemma [S] 
below. Hence all the potentials u'-^\z) axe equal, and all fij = Au(^^/27r are equal 
as well. Finally we get 

Sk'\z) 



Ciz) 



um n^to'- 



lim 



V{z) 



{n - j)Si^\z) ri...(n -k + l)Sn(z) Qk{z) 

for almost all z. This completes the proof of the existence of a compactly supported 
probability measure whose Cauchy transform satisfies (|1.4p in the case when V{z) 
is the limit of a sequence of normalized Van Vleck polynomials of a Heine-Stieltjes 
problem with the leading coefficient of the operator equal to Qk{z). The uniqueness 
of such a measure was obtained in § 

Lemma 9 (see Lemma 8 of [5 ). Let be a sequence of polynomials, such 

that Um '■= degpm — > oo and there exists a compact set K containing the zeros of 
M Pm{z) simultaneously. Finally, let fim md fi'^^ be the root-counting measures of 
Prn{z) and p'^{z) resp. If fim ~> M '^^'^ f-'m ~^ f^' with compact support and u{z) 
and u'{z) are the logarithmic potentials of fi and fi' , then u'{z) < u{z) in the whole 
C. Moreover, u{z) = u'{z) in the unbounded component o/ C \ supp fi. 

Proof. Assume wlog that Pm{z) are monic. Let K he a, compact set containing the 
zeros of all Pm{z). We have 



u{z) 



lim log|p™(z)| 



and 



= lim 



1 



m^oo n^n — 1 

with convergence in L\^^. Hence 

1 



log 



= lim log 



P'ra{z) 



P'miz) 



nmPm{z) 



lim 

m — ^oo 



— = lim log 

m — ^oo ri 

Now, if (f){z) is a positive test function it follows that 
(l){z){u' {z) - u{z))d\{z 



lim log 



1 0(2) log 






J z-C 



< lim 

771 — >00 



(2) 



CI 



■dX{z) < lim 



(j){z)dX{z) 



z-C 
dX{z) < 



where A denotes Lebesgue measure in the complex plane. Since is locally inte- 
grable, the function J (j){z)\z — C,\~^dX{z) is continuous, and hence bounded by a 
constant M for all z in K. Since supp fj,m € K, the last expression in the above 
inequality is bounded by M/um, hence the limit when to — s- oo equals to 0. This 
proves u'{z) < u{z). 

In the complement to supp /i, u{z) is harmonic and u'{z) is subharmonic, hence 
u'{z) — u(z) is a negative subharmonic function. Moreover, in the complement of 
K , p'„{z) I {nmPm{z)) converges uniformly on compact sets to the Cauchy transform 
Cfj,{z) of II. Since Cp_{z) is a non-constant holomorphic function in the unbounded 
component of C \ K, then by the above u'{z) — u{z) — there. By the maximum 
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principle for subharmonic functions it follows that u'{z) — u{z) = holds in the 
unbounded component of C \ supp /i as well. □ 

To accomplish the proof of Theorem [3] we need to show that we have the con- 
vergence for the whole sequence and not just for some subsequence. Assume now 
that the sequence /2„ does not converge to /2. Then we can find a subsequence NN' 
such that fin stay away from some fixed neighborhood of /2 in the weak topology, 
for all n G NN' . Again by compactness, we can find a subsequence NN* of NN' 
such that all the limits for root measures for derivatives exist for j — 0, fc. But 
then p,^^'^ must coincide with p. by the uniqueness and the latter lemma. We get 
a contradiction to the assumption that /i„ stay away from fi for all n G NN' and 
hence all n G NN. □ 



4. Final Remarks 

As an observant reader easily notices the present paper leads to many more 
questions than it provides answers to, some of those already mentioned in the 
introduction. We use this circumstance as an excuse for having this lengthy final 
section. 

I. Call the differential monomial Qk{z)-^ the leading term of a differential operator 

''(-^) = Ei=i Qii^)jli- ^^^^ conjecture states that asymptotically Van Vleck 
and Stieltjes polynomials of the equation (|1.3p behave very similarly to that of the 
equation 

Qk{z)^-^ + V{z)S{z)^0. (4.1) 

Let Polr denote the space of all monic polynomials of degree r. Take a non- 
degenerate c)(z). For some Van Vleck polynomial V{z) of denote by V{z) its 
monic scalar multiple. For a given positive integer n denote by {Vn,i{z)} the set 
of all Van Vleck polynomials V(z) whose Stieltjes polynomials have degree exactly 
n. (Each V{z) is repeated as many times as its multiplicity prescribes, see 
Notice that for sufficiently large n the set {Vn,i{z)} belongs to Polr, i-e. each Van 
Vleck has degree exactly r. Transform now the set {Vn,i{z)} into a finite measure 
CT„(0(z)) in Polr by assigning to each element the finite mass equal to the inverse 
of the cardinality of {Vn^i{z)}. 

Conjecture 1. For any non- degenerate d{z) 

• the sequence cr„(0(z)) converges weakly to a measure I](t)(z)) compactly 
supported in Polr,' 

• the measure T,{d{z)) depends only on the leading monomial Qk{z)-^ of 
5(z). 

Remark 7. One can defined the natural 'projection' of the finite set of polynomials 
cr„(i)(2;)) to the union of the zero loci 3(o'„(c)(2;))) of these polynomials and then 
turn the latter set into a finite measure in the same way. Conjecture [1] implies 
that this set of finite measures converges to the standard measure supported inside 
ConvQ^ , see two left pictures on Fig. U) 

II. Let us present a quite surprising (at least to us) observation about the as- 
ymptotic distribution of Van Vleck versus Stieltjes polynomials obtained through 
numerical experiments. We state it first in the simplest case r = 1. Notice that for 
any nondegenerate Lame operator d{z) — Qk{z) ^^k-i + ■•■ of order k with r — 1 
and sufficently large n there exist exactly (n-|- 1) Van Vleck polynomials all having 
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degree 1, see Proposition [TJ Denote by ^„ the root-counting measure for the union 
of all zeros of these polynomials, and define ^ — lim„_+oo (if it exists). Assuming 
that ^ exists which is strongly supported by our numerics we formulate the following 
fact illustrated on Figure 4 below. 

Observation 1. The support of ^ has a very surprising topological and geometric 
resemblance with that for the limiting root counting measure of the sequence of 
eigenpolynomials for the exactly solvable (i.e. r = 0) operator T = Qk{z)j^, see 
Theorem 3 of jSj. (Notice that in case r = eigenpolynomials are natural analogs 
of Heine-Stieltjes polynomials for such an operator.) 

2r ' ' ' ' ' 2r ' ' ' ' ' 7^ 

1 •■ . 1 • . 

• ■ ' 

-2 '■. -2 

-3L , , , ji , J -3L . . , . ^ 

0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



,4 

Figure 4. The zeros of the eigenpolynomial of degree 40 for the operator Qiz)^^ 
(left) and combined zeros of the 40 linear Van Vlecks having a Stieltjes 
polynomial of degree 39 for the operator j^{Q{z)S{z)) + V{z)S{z) = 
with Q{z) = {z^ + l){z -31- 2){z + 21 - 3), comp. FiglH 

This resemblance seems to persist for larger r. Namely, consider two higher 
Lame operators Di and 02 of the form 



= Qk{z)—i + ■ • ■ and ^2 = Qk{z) 



dz^ dz^ 1 ' 

where AegQk{z) > k and I < k. Let Vn{z) denote the product of all Van Vleck 
polynomials having Stieltjes polynomials of degree n for 1)2 and let 5„(z) denote 
the product of all Stieltjes polynomials of degree n for Oi. Let i/i be the asymptotic 
root-counting measure for the sequence {S'„(z)} and i>2 be the asymptotic root- 
counting measure for the sequence {y„(z)} (if they exist). 

Observation 2. In the above notation the supports of 1^1 and 1^2 have surprising 
geometric and topological similarities, see Figure 5. 

Problem 4. Find an explanation for the above observations. 

III. Curves similar to that in the support of y appeared earlier as Stokes lines 
of linear differential equations with polynomial coefficients. 

Problem 5. What is the exact relation between the WKB-theory for the equation 

y'"' + §4^ = (4.2) 
Q{z) 

and the asymptotic root counting measure /ig^^-j yof the sequence of Stieltjes poly- 
nomials with the sequence of normalized Van Vleck polynomials converging to V ? 

In the case of a generalized Lame equation (ll.ip one at least has a strict definition 
of the global Stokes line for such an equation, see e.g. [7j and the support of the 
corresponding /^^(x) y coincides with a part of this global Stokes line. Even in this 
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0.5 1 1.5 2 2.5 



0.5 1 1.5 2 2.5 



0.5 1 1.5 2 2.5 3 




0.5 1 1.5 2 2.5 3 



Figure 5. Combined zeros of the Stieltjes polynomials for Qiz)j^ and Q{z)jp 
(two left pics) and combined zeros of the Van Vleck polynomials 
for Q(z)^ and Q(2)^ with Q(z) = {z^ + l){z - 31 - 2){z + 21 - 3) 
as described in Observation 2 above. Notice the resemblance of shapes! 



classical case there is an interesting problem of defining what part of the global 
Stokes line is covered by this support. 

For the case of equations of order exceeding 2 the situation is much worse since 
even a good definition of the global Stokes line creates serious problems, see e.g. 
[2], [3], references therein and further publications of the same authors. 

IV. Next question generalizes Problem 3 from the introduction. 

Problem 6. For what plane algebraic eurves in ivith the given coordinate sys- 
tem (z, w) there exists a compactly supported probability measure whose Cauchy 
transform is a section of this algebraic curve almost everywhere in C? 

At the moment the authors know of several such classes of algebraic curves re- 
lated to eigenpolynomials of linear ode, but it is completely unclear how to describe 
the whole set of such curves. 

V. In [3] T. Bergkvist obtained a number of interesting results and conjectures in the 
case of degenerate exactly solvable operators, i.e deg Qk < k + r, see Introduction. 
Motivated by her results we formulate the following conjecture and a question. 

Conjecture 2. For any degenerate Lame operator and any positive integer Nq the 
union of all the roots to polynomials V and S taken over deg S > Nq is always 
unbounded. Therefore, this property is a key distinction between non- degenerate 
and degenerate Lame operators. 

Problem 7. Extend the results of this paper to the case of degenerate higher Lame 
operators. 

VI. Numerical experiments suggest that Stieltjes polynomials of consecutive de- 
grees show a stable root-interlacing pattern along the curves in 6^ where n is the 
corresponding asymptotic root counting measure. This phenomenon is especially 
easy to explain and illustrate in the case of exactly solvable operators, i.e. r = 0. 

Conjecture 3. For any exactly solvable operator d{z) the family {Sn{z)}neN,n>no 
of its eigenpolynomials (deg Sn{z) — n forn CzN) has the interlacing property along 
the support of its asymptotic root counting measure. In other words, the zeros of 
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any two consecutive polynomials Sn+i{z) and Sn{z) interlace along every curve 
segment in 6^ for all sufficiently high degrees n. 

Some caution is required when defining the notion of interlacing since a) the zeros 
of Sn{z) do not he exactly on 6^ and b) 6^ has a nontrivial topological structure. 
One has to remove sufSciently small neighborhoods of the singular points of 6^ 
where several smooth branches meet and to project the zeros of Sn{z) onto 
along some fixed in advance normal bundle to the smooth part of 6^. 

An example illustrating the interlacing phenomenon conjectured above is shown 
on Figure 6. 



1^ 



0.2 



Figure 6. Interlacing pattern for the roots of 6*40(2) and S4i{z) which are the 
eigenpolynomials of the operator z{z — l){z — 

Similar interlacing property was observed for any family of Stieltjes polynomi- 
als of increasing degree as soon as the sequence of their normalized Van Vleck 
polynomials has a limit. Some results about interlacing can be found in [1]. 

VII. In connection with the classical Bochner-Krall problem which asks to describe 
all families of orthogonal polynomials appearing as the families of eigenpolynomials 
for linear differential operators with polynomial coefficients it is natural to ask the 
following. 

Problem 8. Describe all Lame operators for which all their Van Vleck and Stieltjes 
polynomials have only real roots. 

VIII. Finally theorem [3] can be reinterpreted as the statement that the level curves 
of the logarithmic potential Ufi{z) of the measure which we construct in this 

theorem are trajectories of the differential ( — 1 ) q^^I~^ dz'' of order k. Since almost 
all level curves of are closed curve we are tempted to call the latter differential 

Strebel. This makes perfect sense when k = 2, see [13] and [H]. But for k > 2 even 
a definition of a Strebel differential is so far missing. 

Problem 9. Define a notion of a Strebel (rational) differential of order k > 2. 

References 

[1] A. Agncw, A. Bourget, T. McMillen, On the zeros of complex Van Vleck polynomials, J. 
Comp. Appl. Math, 223, (2009), 862-871. 

[2] T. Aoki, T. Kawai and Y. Takei, New turning points in the exact WKB analysis for higher 
order ordinary diffrential equations, Analyse algebrique des perturbations singulicrcs. I, Her- 
mann, (1994), 69-84. 



ON HIGHER HEINE-STIELTJES POLYNOMIALS 



17 



T. Aoki, T. Kawai and Y. Takei, On the exact WKB analysis for the third order ordinary 
differential equations with a large parameter, Asian J. Math., vol 2, (1998), 625—640. 
T. Bergkvist, On asymptotics of polynomial eigenf unctions for exactly solvable differential 
operators, J. Approx. Theory 149(2), (2007), 151-187. 

T. Bergkvist and H. RuUgard, On polynomial eigenfunctions for a class of differential oper- 
ators, Math. Res. Lett., voh9, (2002), 153-171. 

B. Shapiro, Algebro-geometric aspects of Heine-Stieltjes theory. rarXiv:0812.4193] 

M. Fedoruyk, Asymptotic analysis. Linear ordinary differential equations. Translated from 

the Russian by Andrew Rodick. Springer- Verlag, Berlin, (1993), viii-|-363 pp. 

E. Heine, Handbuch der Kugelfunctionen, vol.1, Berlin: G. Reimer Verlag, (1878), 472-479. 

T. Hoist, B. Kruglikov, B. Shapiro, Local study of multi-valued analytic differentials and their 

trajectories, in preparation. 

J. A. Jenkins, Univalent functions and conformal mapping. Ergebnisse der Mathematik und 
ihrer Grenzgebiete. Neue Folge, Heft 18. Reihe: Modernc Funktionentheorie Springer- Verlag, 
Berlin-Gttingen-Heidelberg, (1958), vi-|-169 pp. 

A. Martfnez-Finkelshtein, E. A. Rakhmanov, Critical measures, quadratic differentials, and 
weak limits of zeros of Stieltjes polynomials, arXiv:0902.0193 

A. Martmez-Finkelshtein and E. SafT, Asymptotic properties of Heine-Stieltjes and Van Vleck 
polynomials, J. Approx. Theory, vol. 118(1), (2002), 131-151. 

G. Polya, Sur un theoreme de Stieltjes, C. R. Acad. Sci Paris vol. 155, (1912), 767-769. 

B. Shapiro, K. Takemura, M. Tater, On spectral polynomials of the Heun equation. II, in 
preparation. 

T. Stieltjes, Sur certains polyndmes qui verifient une equation differentielle lineaire du second 

ordre et sur la theorie des fonctions de Lame, Acta Math. vol. 8, (1885), 321-326. 

K. Strebel, Quadratic differentials. Ergebnisse der Mathematik und ihrer Grenzgebiete (3) 

[Results in Mathematics and Related Areas (3)], 5. Springer- Verlag, Berlin, (1984), xii-|-184 

pp. 

K. Strebel, On the metric |/(^)|'''|d^| with holomorphic f. Complex analysis (Warsaw, 1979), 
Banach Center Publ., 11, PWN, Warsaw, (1983), 323-336. 

G.Szego, Orthogonal polynomials. Fourth edition. American Mathematical Society, Collo- 
quium Publications, Vol. XXHL American Mathematical Society, Providence, R.L, 1975. 
xiii+432 pp. 

E. Van Vleck, On the polynomials of Stieltjes, Bull. Amer. Math. Soc. vol. 4, (1898), 426-438. 
E. Whittaker and G. Watson, A course of modern analysis. An introduction to the general 
theory of infinite processes and of analytic functions; with an account of the principal tran- 
scendental functions. Reprint of the fourth (1927) edition. Cambridge Mathematical Library. 
Cambridge University Press, Cambridge, 1996. vi-|-608 pp. 

Department of Mathematics, Stockholm University, SE-106 91 Stockholm, Sweden 
E-mail address: holstOmath. su. se 



Department of Mathematics, Stockholm University, SE-106 91 Stockholm, Sweden 
E-mail address: shapiroOmath. su. se 



